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Abstract. The shell model Monte Carlo (SMMC) approach provides a powerful method for 
the microscopic calculation of statistical and collective nuclear properties in model spaces that 
are many orders of magnitude larger than those that can be treated by conventional methods. 
We discuss recent applications of the method to describe the emergence of collectivity in the 
framework of the configuration-interaction shell model and the crossover from vibrational to 
rotational collectivity in families of rare-earth nuclei. We have calculated state densities of these 
rare-earth nuclei and find their collective enhancement factors to be correlated with the pairing 
and shape phase transitions. We also discuss an accurate method to calculate the ground-state 
energy of odd-even and odd-odd nuclei, circumventing the sign problem that originates in the 
projection on an odd number of particles. We have applied this method to calculate pairing 
gaps in families of isotopes in the iron region. 



1. Introduction 

A major challenge in nuclear theory is the derivation of nuclear properties from the underlying 
effective nuclear interactions. Advances in many-body methods and high-performance 
computing have enabled significant progress in the application of ab initio methods such as 
Green's function Monte Carlo [1], and the no-core shell model [2] and its symplectic extension [3]. 
However, these methods are limited to light nuclei. Density functional theory (DFT) has made 
much progress and is unique in providing a global theory of nuclei [1] but correlations beyond 
mean-field theory are often important and require extensions of DFT [5]. 

The configuration-interaction (CI) shell model approach provides an attractive alternative as 
it accounts for both shell effects and correlations [6] . With the help of increasingly more efficient 
matrix diagonalization algorithms, the CI approach has been successful in describing the low- 
energy spectroscopy of /p-shell nuclei. However, in heavier nuclei conventional diagonalization 
methods become limited by the combinatorial growth of the dimensionality of the many-particle 
space, hindering applications of the conventional CI approach in medium-mass and heavy nuclei. 

The shell model Monte Carlo (SMMC) method [5] has enabled calculations in model 
spaces that are many orders of magnitude larger than those that can be treated by conventional 
methods. Monte Carlo methods for fermions are often limited by the so-called sign problem. 



However, the dominant components [9] of realistic effective nuclear interactions have a good sign 
and are often sufficient for realistic calculations of statistical and collective properties of nuclei. 
The smaller, bad sign components can be treated by the method of Ref. |8j. 

Here we discuss recent developments in SMMC, and in particular the application of SMMC 
in the study of collective and statistical properties of heavy nuclei \10\ lllj . Another recent 
progress has been the introduction of a method to calculate the ground-state energy of odd-even 
and odd-odd nuclei, circumventing a sign problem that originates from the projection on an odd 
number of particles [12j . 



2. The shell model Monte Carlo approach 

2.1. Hubbard- Stratonovich transformation 

The SMMC is based on the Hubbard-Stratonovich (HS) transformation |13] . in which the Gibbs 
operator e~@ H of a system with an Hamiltonian H at an inverse temperature (3 = 1/T is 
expressed as a functional integral 

e-t> H = [ D[a]G a U a (1) 



over external auxiliary fields <t(t). Here G a is a Gaussian weight and U a is a one-body propagator 
describing non-interacting nucleons moving in external fields <t(t). 

Using the HS transformation, the thermal expectation value of an observable O can be written 

as 

Tr(Oe-^) _ jD[a]WMO) a 
{ > Tr(e-f> H ) fD[a}W a <f> a ' { ' 

where W a = Gr CT |Tr U a \ is a positive-definite function and <J? CT = Tr U a /\Ti U a \ is known as 
the Monte Carlo sign function. The quantity {0) a = Tr (OU&) /Tr U a describes the thermal 
expectation value of the observable at a given configuration of the auxiliary fields a. 

In the SMMC approach H is a rotationally invariant shell-model Hamiltonian defined in a 
truncated many-particle space spanned by a specified set of N s single-particle orbitals (this set 
can be different for protons and for neutrons). The quantities in the integrands of Eq. ([2]) can 
be expressed in terms of the N s x N s matrix JJ a representing U a in the single-particle space. 
For example, the grand-canonical trace of U a is given by Tr U a = det(l + U CT ). 

In the finite nucleus it is important to use the canonical ensemble with fixed number of 
protons and neutrons. This is done using particle-number projection, which for a finite number 
of orbitals can be expressed as a discrete Fourier sum. For example, the trace of U a at fixed 
particle number A is given by 

p —/3fJ.A N s 

Tr A U a = -j^- e ~ l4> det (l + e^™e^U CT J , (3) 

s m=l 

where A is the particle-number operator, <j) m = 2nm/N s (m = 1, . . . , N s ) are quadrature points 
and \x is a real chemical potential [required for the numerical stabilization of the sum in (J3J)]. 

2.2. Monte Carlo method 

The imaginary-time interval (0, (3) is discretized in steps of Af3. The integration over the 
auxiliary fields is of very high dimension and is carried out by a Monte Carlo method. The 
integration over each a field at a given time is approximated by a three-point quadrature formula 
and importance sampling is used to choose M uncorrelated configurations according to the 
positive-definite weight function W a . The expectation value in ([2]) is then estimated from 

(O) » Hk(Q)° k ®<Tk / 4 ) 



3. Collectivity in the configuration-interaction shell model approach 

Heavy nuclei are known to exhibit collectivity in their low-energy spectra and electromagnetic 
transition rates. Phenomenological nuclear structure models describe well the various types 
of collectivity such as vibrational collectivity in spherical nuclei and rotational collectivity 
in deformed nuclei. However, much less is understood about the microscopic emergence of 
collectivity in the framework of a CI shell model approach. A shell model treatment of heavy 
deformed nuclei requires very large model spaces and that are accessible by the SMMC method. 

We have extended the SMMC method to heavy nuclei using a proton-neutron formalism in 
which the protons and neutrons can occupy different shells [TU]. In particular, we have applied 
the method to rare-earth nuclei using the 50 — 82 major shell plus the I/7/2 orbital for protons 
and the 80 — 126 major shell plus the 0h u / 2 and lgg/2 orbitals for neutrons. 

The SMMC method enables us to calculate thermal and ground-state observables in very 
large model spaces, but it does not provide detailed spectroscopic information. Since nuclear 
collectivity is often identified by its spectroscopic characteristics, it is necessary to find a thermal 
observable that is sensitive to the specific type of collectivity. We have identified such an 
observable in (J 2 )t ; where J is the total angular momentum of the nucleus [10\ 111] . For an 
even-even nucleus and at sufficiently low temperatures, this observable is given by 

-E 2+ /T 

30- r-r vibrational band 



<J>«< (5) 
■^r—T rotational band 

where E 2 + is the excitation energy of the first 2 + level. 



3.1. Rotational collectivity in heavy deformed nuclei 




Figure 1. (a) (J 2 ) T and (b) thermal energy E(T) versus T for 162 Dy. The SMMC results 
(solid circles with error bars) are compared in (a) with fits to the rotational model (solid line) 
and vibrational model (dotted-dashed line) [see Eq. (jSJ], and in (b) with the rotational model 
(solid line). The dashed lines are obtained from the lowest five experimental rotational bands. 
Adapted from Ref. [ID] . 



An important issue is whether we can describe the rotational collectivity of heavy deformed 
nuclei in the framework of a truncated spherical shell model space. In Ref. [TO] we studied 
a strongly deformed rare-earth nucleus 162 Dy using the SMMC method and demonstrated its 
rotational character through the low-temperature behavior of (J 2 )t and the thermal energy 
E(T). Figure [I] shows these observables as a function of temperature. For comparison we 
also show for (J 2 )t the fits to the rotational (solid line) and vibrational (dotted-dashed line) 
models in Eq. ([5]). The SMMC results clearly indicate the rotational character of 162 Dy in 



the truncated shell model space. Furthermore, the fit to the rotational model gives a ground- 
state moment of inertia of / g . s . = 35.5 ± 3.3 MeV -1 , in agreement with the experimental value 



of /, 



g.s. 



3/E 2 + = 37.2 MeV -1 . The solid line in panel (b) is a fit to the rotational model 



prediction for the thermal energy E(T) = Eq + T, where Eq is the ground-state energy. 
3.2. Crossover from vibrational to rotational collectivity 
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Figure 2. (J 2 ) T function of temperature in the even 148 154 Sm isotopes. The SMMC results 
(open circles with error bars) are compared with Eq. (|6|) using a complete set of experimentally 
known discrete nuclear levels and an experimentally determined BBF state density (except for 
154 Sm, for which only discrete nuclear levels are used). Adapted from Ref. [TT] , 



Even samarium and neodymium isotopes display a shape phase transition from spherical to 
deformed nuclei as the number of neutrons increases from shell closure towards the mid-shell 
region. In the SMMC approach, a signature of this transition can be observed in the low- 
temperature behavior of (J 2 )t- In Fig. [2]we compare in the even 148_154 Sm isotopes the SMMC 
results (circles) for (J 2 )t with results that are derived from the experimental data (solid lines) 
according to 

1 ( ' SL. r°° \ 

(J 2 )t = (£ Ji(Ji + 1)(2J* + l)e- E */ T + a!E x p(E x ) <J 2 }^ e~ E */ T ) . (6) 

Here Z(T) = Y.i { 2 Ji + l)e~ El/T + Je dE x p(E x )e- Ex/T is the corresponding experimental 
partition function. The summations in %]) and in the expression for Z{T) are carried over a 
complete set of experimentally known low- lying levels i with excitation energy Ei and spin Jj, 
and p(E x ) is the back-shifted Bethe formula (BBF) with parameters that are determined from 
level counting and neutron resonance data. (J 2 )e x (the average value of J 2 at a given excitation 
energy E x ) is calculated from (J 2 )^ = 3<r 2 , where a 2 is the rigid-body spin-cutoff parameter. 
We find good agreement between the SMMC results and the experimental estimates. 148 Sm is 
spherical in its ground state and its (J 2 )t exhibits a "soft" response to temperature, typical 
of vibrational nuclei. In contrast, (J 2 )t in 154 Sm exhibits a linear response to temperature, 
indicating the rotational character of this nucleus. A similar study of even 144_152 Nd isotopes 
confirms that 144 Nd and 146 Nd are vibrational, while 150 Nd and 152 Nd are rotational. 



4. State densities in heavy nuclei 

4-1. State densities 

The SMMC method has been successful in the calculation of state densities in medium-mass 
nuclei |14tll5|. In the SMMC method we calculate the average thermal energy E(f3) as a function 
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Figure 3. Total state densities of the even 148_154 Sm isotopes. The SMMC results (open circles) 
are compared with the level counting data (histograms), the neutron resonance data (triangles), 
the experimentally parametrized BBF state densities (solid lines) and the HFB densities (dashed 
lines). Arrows denote proton and neutron pairing phase transitions and the thick arrow denotes 
a shape transition. Adapted from Ref. 



of /3, from which we can compute the canonical partition function. The average state density is 
then determined by applying the saddle-point approximation to the integral that expresses the 
state density as the inverse Laplace transform of the canonical partition function. 

In Fig. [3] we show the SMMC state densities (open circles) in the even 148_154 Sm isotopes 
and compare them with level counting data at low energies (histograms) and neutron resonance 
data (triangles) when available. The solid lines are BBF state densities that are determined from 
fits to experimental data. The dashed lines are the densities obtained in the finite-temperature 
Hartree-Fock-Bogoliubov (HFB) approximation. The HFB results describe the intrinsic density 
of states, and the enhancement observed in the SMMC state density arises from collective 
states (vibrational and rotational) that are built on top of the intrinsic states. In the HFB 
approximation we observe proton and neutron pairing phase transitions (arrows) and shape 
phase transitions (thick arrows). 148 Sm is spherical in its ground state and exhibits pairing 
transitions only, while 150 Sm is deformed in its ground state and thus undergoes also a shape 
transition to a spherical shape at E x ~ 12.5 MeV. Shape phase transitions also occur in 152 Sm 
and 154 Sm at excitation energies that increase with mass number (not shown in the figure). 



4-2. Collective enhancement factors 

Collective effects in the nuclear state density are difficult to model microscopically and are 
usually accounted for by phenomenological vibrational and rotational collective enhancement 
factors [TS]. In Ref. [TT] we have proposed to use the ratio of the SMMC and the HFB state 
densities, K = psmmc/phfb, as a microscopic measure of the collective enhancement factor. In 
Fig. H] we show K as a function of the excitation energy E x for the same samarium isotopes 
shown in Fig. [3J In 148 Sm, which is spherical in its ground state, any collectivity should be 
exclusively due to vibrational excitations, whereas the heavier deformed samarium nuclei are 
expected to have both vibrational and rotational collectivity. We observe that all collectivity 
disappears in 1 8 Sm (i.e., K ss 1) at excitation energies above the pairing transitions, suggesting 
that the decay of the vibrational collectivity is associated with the pairing phase transitions. In 
the heavier samarium isotopes K reaches a local minimum above the pairing transition energies 
but increases at higher energies until it reaches a plateau and then decays to unity in the vicinity 
of the shape transition energy. This behavior can be explained by the contribution of rotational 
collectivity to K that decays to unity once the nucleus becomes spherical and can no longer 
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Figure 4. Total collective enhancement factor K (see text) in the even 148 154 Sm isotopes as 
a function of excitation energy E x . Arrows are as in Fig. [3l Adapted from Ref. [llj . 



support rotational bands. 



5. Odd-even and odd-odd nuclei 

The projection on an odd number of particles leads to a sign problem in SMMC (even for good 
sign interactions). This problem is particularly severe at low temperatures and makes it difficult 
to calculate the ground-state energy of odd-even and odd-odd nuclei. We have developed a 
method that circumvents this sign problem for the ground-state energy by using asymptotic 
properties of the single-particle Green's functions of the neighboring even-even nuclei [12] . 



5.1. Green's function method 

For a rotationally invariant Hamiltonian, the scalar imaginary-time single-particle Green's 
functions are given by 

TY.A e-^TE (o) 
Gu{r) = l - „„ ] - , (7) 



TY4 e~P R 

where v = (nlj) labels the nucleon single-particle orbital with radial quantum number n, orbital 
angular momentum I and total spin j. Here T denotes time ordering and a^ m (r) = e rH a vm e~ 
is an annihilation operator of a nucleon at imaginary time r (— ft < r < j3) in a single-particle 
state with orbital v and magnetic quantum number m. 

Consider an even-even nucleus A = (Z, N), with Z protons and N neutrons. Assuming that 
the ground state of this nucleus has spin 0, the r-dependence of the Green's function has the 
asymptotic form G v {t) ~ e _A£ ' J=j( '" 4± ^' r ', where A± denote the even-odd nuclei (Z,N±1) when 
v is a neutron orbital and the odd-even nuclei [Z ± 1,N) when v is a proton orbital. The + 
and — subscripts should be used for r > and r < 0, respectively. AEj = j(A±) is the difference 
between the energies of the lowest spin J eigenstate of the .4±-particle nucleus and the ground 
state of the ^.-particle nucleus. In this asymptotic regime for r, we can calculate AEj(A±) from 
the slope of In G v (t). The minimum of AEj(A±) over all possible values of j gives the difference 
between the ground-state energy of the A± nuclei, E gs (A±), and the ground-state energy of the 
A nucleus, E gs (A). Both E gs (A) and G u (t) are related to the even-even nucleus and can be 
calculated in SMMC without a sign problem (for a good sign interaction). 
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Figure 5. Neutron pairing gaps A n as a function of mass number A in /p^g^-shell nuclei. 
The gaps calculated with the Green's function method (solid circles) are compared with the 
experimental gaps (open circles). Adapted from Ref. [12] . 



5.2. Pairing gaps in iron-region nuclei 

The method of Sec. 15. H is particularly useful for calculating pairing gaps (i.e., odd-even staggering 
of masses) because the odd-even ground state energy differences are calculated directly. In 
Fig. [5]we compare the SMMC pairing gaps of iron-region nuclei with their experimental values. 
We have used the complete /pgg/2-shell model space with the Hamiltonian of Ref. |14j . The 
statistical errors for the SMMC pairing gaps (~ 0.01 MeV) are not visible in the figure. The 
theoretical gaps are in reasonable agreement with the experimental values. 



6. Conclusion 

We have discussed recent progress in the SMMC approach. We have shown that the crossover 
from vibrational to rotational collectivity in rare-earth nuclei can be described microscopically 
in the framework of the CI shell model approach. We have also described a method to calculate 
the ground-state energy of odd-even and odd-odd nuclei that circumvents the odd-particle sign 
problem and applied it to calculate pairing gaps of nuclei in the iron region. 
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